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Aromatic  cations  such  as  imidazolium  and  its  derivatives  have 
delocalized  charges  and,  therefore,  are  often  used  to  synthesize 
functional  nanomaterials,  low-melting  salts  and  ionic  liquids\ 
Recently,  energetic  low-melting  salts  and  ionic  liquids  based  on 
triazolium  and  its  derivatives  have  received  considerable 
attention^'"^.  A  characteristic  property  of  aromatic  compounds  is 
their  ability  to  form  n-n  stacking  structures,  which  can  have 
significant  influence  on  properties.  Cation-cation  n-n  stacking 
structures  are  less  common,  but  have  been  observed  in 
organometallic  crystals^  and  imidazolium  salts^.  Cation-cation  n- 
n  stacking  dimers  can  be  formed  in  a  dielectric  solvent  due  to 
significant  screening  of  the  charge-charge  repulsion  forces.^ 
However,  such  species  in  small  clusters  have  not  been  reported. 
Calculations  reported  in  this  work  suggest  that  cation-cation  n-n 
stacking  structures  can  be  formed  between  two  1,2,4-triazoiium 
cations  in  small  clusters  of  the  two  cations  and  two  anions. 

The  ab  initio  calculations  were  performed  with  the  electronic 
structure  package  GAMESS  Geometries  were  optimized  using 
second  order  perturbation  theory  (MP2)  with  the  aug-cc- 
pVDZ^^  basis  set.  Single  point  energies  were  calculated  for 
selected  structures  using  MP2  with  the  cc-pVDZ,  aug-cc-pVDZ, 
aug-cc-pVTZ  and  aug-cc-pVQZ  basis  sets^^,  and  coupled  cluster 
singles,  doubles  with  perturbative  triples  [CCSD(T)]  with  the 
cc-pVDZ  and  aug-cc-pVDZ  basis  sets. 
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Figure  1.  MP2/aug-cc-pVDZ  optimized  structures  of  two  pairs  of  1,2,4- 
triazoiium  (1,2,4-triazole)  and  dinitramide  (dinitramine)  molecules.  H  is 
white,  C  is  gray,  O  is  red,  N  is  blue. 


Numerous  structures  can  be  found  for  two  pairs  of  1,2,4- 
triazoiium  and  dinitramide  ions  or  two  pairs  of  1,2,4-triazole  and 
dinitramine  neutral  molecules.  An  extensive  geometry  search  was 
performed  for  the  ionic  and  neutral  structures.  The  lowest  energy 
MP2/aug-cc-pVDZ  structures  are  shown  in  Figure  1.  In  the  most 
stable  ionic  two-pair  structure,  each  1,2,4-triazoiium  forms  two 
hydrogen  bonds,  via  the  hydrogens  on  the  N  atoms,  to  the  O 
atoms  of  the  dinitramide  ions.  Interestingly,  this  structure  exhibits 
n-n  stacking  of  the  two  cationic  1,2,4-triazoiium  rings.  The 
interplane  distance  is  ~3.2  A,  with  a  parallel  displacement  of -1.4 
A.  The  nitro  groups  of  the  dinitramide  ions  are  approximately 
perpendicular  to  each  other.  In  the  most  stable  neutral  two-pair 
structure,  each  1,2,4-triazole  forms  one  N-H-N  hydrogen  bond  to 
a  dinitramine.  There  is  n-n  stacking  between  the  two 
approximately  parallel  1,2,4-triazole  rings,  with  an  interplane 
distance  of  -3.2  A  and  a  parallel  displacement  of  -1.4  A.  An 
ionic  structure  was  created  from  the  neutral  structure  by 
repositioning  protons  onto  the  1,2,4-triazole  molecules.  Upon 
geometry  optimization,  the  protons  transfer  back  to  the 
dinitramide  ions  to  form  neutral  pairs  with  no  energy  barrier. 

It  is  of  interest  to  determine  the  cluster  size  at  which  the  ion 
pair  structures  become  more  stable  than  the  corresponding  neutral 
pair  structures.  A  previous  study  predicted  that  ion  pair  dimers  are 
typically  higher  in  energy  than  neutral  pair  dimers."^  The 
MP2/aug-cc-pVDZ  energy  of  the  ionic  tetramer  in  Figure  1  is  1.5 
kcal/mol  lower  than  that  of  the  neutral  one.  Including  zero  point 
energies  (ZPE)  obtained  with  MP2/aug-cc-pVDZ  harmonic 
vibrational  frequencies,  the  ionic  tetramer  is  1 .2  kcal/mol  lower  in 
energy  than  that  of  the  neutral  one.  The  MP2  method  tends  to 
predict  higher  energies  for  ionic  species  vs.  neutral  species."^ 
Since  the  computational  cost  of  CCSD(T)/aug-cc-pVDZ  energy 
calculations  of  these  two  tetramer  structures  is  prohibitive,  these 
energies  were  approximated  from  the  MP2/aug-cc-pVDZ  energies 
by  estimating  the  electron  correlation  energy  differences  using 
three  independent  methods:  (1)  the  differences  between  the 
MP2/cc-pVDZ  and  CCSD(T)/cc-pVDZ  energies  of  the  tetramers, 
(2)  the  differences  between  the  MP2/aug-cc-pVDZ  and 
CCSD(T)/aug-cc-pVDZ  energies  of  the  twelve  pairs  of  dimers  in 
these  two  tetramers,  and  (3)  the  differences  between  the 
MP2/aug-cc-pVDZ  and  CCSD(T)/aug-cc-pVDZ  energies  of  the 
eight  monomers  in  these  two  tetramers.  Using  these  three 
methods,  the  estimated  CCSD(T)/aug-cc-pVDZ  energy  of  the 
ionic  tetramer  is  lower  than  that  of  the  neutral  tetramer  by  6.0,  7.6 
and  8.0  kcal/mol,  respectively.  Including  ZPE,  these  three  values 
become  5.7,  7.3  and  7.7  kcal/mol,  respectively.  Further 
corrections  for  the  basis  set  effects  from  aug-cc-pVDZ  to  aug-cc- 
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pVTZ  and  to  aug-cc-pVQZ  are  small.  One  ean  therefore  eonelude 
that  for  two  pairs  of  1,2,4-triazolium  and  dinitramide,  ionie 
struetures  are  more  stable  than  the  eorresponding  neutral 
struetures. 

The  driving  foree  for  the  formation  of  the  unusual  eation-eation 
71-71  staeking  strueture  is  also  very  interesting.  Geometry 
optimization  of  the  ionie  strueture  using  HF,  whieh  does  not 
eontain  eleetron  eorrelation,  results  in  an  open-book  strueture  in 
whieh  the  two  1,2,4-triazolium  planes  form  a  dihedral  angle  of 
-120°  (no  n-n  staeking).  Clearly,  eleetron  eorrelation  (a  major 
part  of  the  aromatie  n-n  attraetion)  in  the  MP2  method  is  the 
driving  foree  for  the  formation  of  the  eation-eation  n-n  staeking 
strueture  from  the  open-book  non-staeking  strueture. 

The  tetramer  interaetion  energy  and  the  six  isolated  dimer 
interaetion  energies  of  the  four  monomers  in  the  ionie  tetramer 
strueture  in  Figure  1  were  derived  by  taking  the  eorresponding 
energy  differenees  between  the  tetramer  and  isolated  dimers  and 
monomers,  with  no  eorreetion  for  basis  set  superposition  errors. 
The  same  geometries  as  in  the  tetramer  were  used  for  the  isolated 
dimers  and  monomers,  and  the  aug-ee-pVDZ  basis  set  was  used. 
Using  HF,  the  tetramer  interaetion  energy  is  -217.1  keal/mol;  the 
repulsion  energies  between  the  two  1,2,4-triazolium  eations  and 
the  two  dinitramide  anions  are  +77.1  and  +52.9  keal/mol, 
respeetively;  the  attraetion  energies  between  the  two  eations  and 
two  anions  are  -96.8,  -89.3,  -89.3,  -96.8  keal/mol  (symmetrie  due 
to  the  symmetry  of  the  tetramer  strueture).  The  -217.1  keal/mol 
tetramer  interaetion  energy  is  roughly  the  sum  (-242.2  keal/mol) 
of  the  dimer  interaetions,  with  a  +25.1  keal/mol  differenee  due  to 
many-body  effeets.  These  HF  interaetion  energies  are  dominated 
by  Coulomb  interaetions.  For  example,  a  Redueed  Variational 
Spaee^^  analysis  of  the  two  1,2,4-triazolium  eations  leads  to  the 
following  terms  (in  keal/mol):  total  =  +77.1,  Coulomb  =  +76.9, 
exehange-repulsion  =  +4.9,  polarization  =  -3.8,  eharge  transfer  =  - 
0.5  and  mixing  terms  -0.4.  In  general,  the  driving  foree  for  the 
formation  of  ionie  elusters  is  the  Coulomb  attraetion  between  the 
opposite-eharge  ions,  whieh  is  typieally  mueh  stronger  than  the 
repulsion  between  the  same-eharge  ions. 

However,  as  shown  above,  additional  dispersion  attraetion  due 
to  eleetron  eorrelation  is  neeessary  for  the  formation  of  the 
eationie  n-n  staeking  strueture  from  the  (HF)  open-book  non- 
staeking  strueture.  The  MP2  tetramer  interaetion  energy  is  -266.7 
keal/mol,  49.6  keal/mol  stronger  than  the  HF  one,  due  to  the 
eleetron  eorrelation  attraetion;  the  MP2  interaetion  is  roughly  the 
sum  (-294.5  keal/mol)  of  the  MP2  dimer  interaetions,  differing  by 
+27.8  keal/mol  due  to  many-body  effeets.  The  CCSD(T)  tetramer 
interaetion  energy  is  not  available;  the  sum  of  the  CCSD(T)  dimer 
interaetions  is  -290.9  keal/mol,  very  elose  to  the  MP2  one.  For  the 
two  1,2,4-triazolium  eations,  the  HF,  MP2  and  CCSD(T) 
interaetion  energies  are  +77.1,  +69.3  and  +70.4  keal/mol, 
respeetively;  the  eleetron  eorrelation  eorreetions  to  the  HF 
interaetion  energy  predieted  by  MP2  and  CCSD(T)  are  -7.8  and  - 
6.7  keal/mol,  respeetively. 

Based  on  the  strueture  of  two  1,2,4-triazolium  -  dinitramide  ion 
pairs  (Figure  1),  an  initial  strueture  of  two  1,2,4-triazolium  - 
ehloride  ion  pairs  was  eonstrueted,  and  then  optimized  using 
MP2/aug-ee-pVDZ.  In  the  optimized  strueture  shown  in  Figure  2, 
eaeh  1,2,4-triazolium  ion  forms  two  N-H-Cl  hydrogen  bonds. 
Apparently,  n-n  staeking  of  the  two  eationie  1,2,4-triazolium 
rings  is  formed  in  this  strueture  in  spite  of  the  small  sizes  of  the 
ehloride  ions.  The  interplane  distanee  is  -2.9  A,  with  a  parallel 
displaeement  of -1.0  A.  The  strueture  is  not  symmetrie  due  to  the 
parallel  up/down  displaeement  (front  view  in  Figure  2).  Clearly, 
very  strong  N-H-Cl  hydrogen  bonds  are  formed  in  this  strueture 
so  that  the  four  H  atoms  bend  inward.  A  neutral  strueture  of  two 


pairs  of  1 ,2,4-triazole  and  hydrogen  ehloride  (HCl)  moleeules  was 
eonstrueted  starting  from  the  1,2, 4- triazole  -  dinitramine  neutral 
tetramer  in  Figure  1,  and  then  optimized  with  MP2/aug-ee-pVDZ. 
The  neutral  strueture  optimized  to  the  ionie  strueture  shown  in 
Figure  2,  with  no  energy  barrier.  Again,  eleetron  eorrelation 
attraetion  is  the  driving  foree  for  the  formation  of  the  n-n  staeking 
strueture:  HF  geometry  optimization  of  the  strueture  leads  to  a 
planar  strueture  in  whieh  all  the  atoms  are  exaetly  in  the  same 
plane. 


Top  view  Front  view 


Figure  2.  MP2/aug-cc-pVDZ  optimized  structures  of  two  pairs  of  1,2,4- 
triazolium  chloride  ions.  H  is  white,  C  is  gray,  N  is  blue.  Cl  is  green. 

In  eonelusion,  ab  initio  ealeulations  suggest  that  eation-eation 
n-n  staeking  struetures  ean  exist  in  very  small  ionie  elusters  sueh 
as  two  1,2,4-triazolium  eations  and  two  dinitramide  or  ehloride 
anions.  The  strueture  motifs  and  interaetion  patterns  found  in  this 
study  provide  new  understanding  of  ionie  materials  with  aromatie 
rings. 
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The  existence  of  cation-cation  n-n  stacking  in  1,2,4-triazoliunn  -  dinitramide  tetramer  and  1,2,4-triazoliunn  -  chloride  tetramer 
(two  cations  and  two  anions)  is  predicted  based  on  the  structures  optimized  using  second  order  perturbation  theory  (MP2).  In 
the  most  stable  tetramer  structure  of  1 ,2,4-triazolium  dinitramide,  n-n  stacking  is  formed  with  an  interplane  distance  of  ~3.2  A 
and  a  parallel  displacement  of  -1.4  A.  In  the  most  stable  tetramer  structure  of  1 ,2,4-triazolium  chloride,  n-n  stacking  is  formed 
with  an  interplane  distance  of  -2.9  A  and  a  parallel  displacement  of  -1 .0  A. 
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